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Abstract 


Early versions of the Langley nuclear fragmentation code NUCFRAG 
(and a publicly released version called HZEFRGl) assumed straight-line 
trajectories throughout the interaction . As a consequence, NUCFRAG 
and HZEFRGl give unrealistic cross sections for large mass removal 
from the projectile and target at low energies. A correction for the 
distortion of the trajectory by the nuclear coulomb fields is used herein 
to derive fragmentation cross sections. A simple energy-loss term is 
applied to estimate the energy downshifts that greatly alter the coulomb 
trajectory at low energy . The results, which are far more realistic 
than prior versions of the code, should provide the data base for future 
transport calculations. The systematic behavior of charge-removal cross 
sections compares favorably with results from low-energy experiments. 


Introduction 

The cross sections in the radiation physics pro- 
gram at the Langley Research Center have a long 
history. Early transport calculations (ref. 1) were de- 
veloped using the work of Waddington and coworkers 
(ref* 2) that was derived mainly from nuclear emul- 
sion data obtained in high-altitude balloons. These 
cross sections were soon abandoned in favor of cor- 
rections to Rudstam’s work (ref. 3) that were derived 
by Silberberg and Tsao (ref. 4) from radionuclide 
production data. The Silberberg-Tsao model was 
mainly appropriate for high charge and energy ions 
(HZE) fragmentation on hydrogen targets, although 
a modification of the scaling procedure of Lindstrom 
et al. (ref. 5) was used for arbitrary target nuclei 
(ref. 6). The cross sections of reference 6 were used 
to explain the Bragg curves of Ne beams in labo- 
ratory experiments (ref. 7). The scaled Silberberg- 
Tsao cross sections (augmented with Bertini (ref. 8) 
cross sections for the lightest fragments n , p , d , t, 
he, and a) significantly underestimated the experi- 
mental Bragg curves. The cause, in part, resulted 
from the lack of mass and charge conservation of the 
cross-section data set. Cross-section renormalization 
to ensure mass and charge conservation improved the 
agreement (ref. 7). 

The next stage of transport development was 
based on a study of atmospheric air shower data 
(refs. 9 and 10). No agreement could be found with 
the charge 3-10 group as measured in the atmosphere 
because of the very large cross sections predicted by 
reference 6 for the production of Li, Be, and B by 
air nuclei from galactic cosmic rays incident at the 
top of the atmosphere (ref. 10). The production of 
a version of the model by Bowman, Swiatecki, and 
Tsang (ref. 11) was attempted with some improve- 
ment, but unsatisfactory 7 agreement was found with 
the air shower data (ref. 10). At that point, the 
first version of the NUCFRAG model was derived 


as a modification of reference 11, and good agree- 
ment was obtained with the air shower data (ref. 10). 
The NUCFRAG code calculates the average excita- 
tion from frictional forces at the interface of the inter- 
action zone and an empirical correction to the surface 
energy for highly misshapened nuclei (ref. 10). 

The frictional forces are derived from two-body 
collisional processes, and the corresponding excita- 
tion energy is a fluctuating variable. The average 
value used in reference 10 was replaced by a random 
variable in reference 12 and was assumed to fluctu- 
ate between its maximum and minimum values with 
equal probability. This procedure gave a satisfac- 
tory fit to the fragmentation data of Westfall et al. 
(ref. 13) and Heckman (ref. 14). 

More recently, the matter densities were used to 
determine nuclear radii, the energy- dependent nu- 
clear mean free path was added, and an improved 
fit of the value of empirical correction to the exci- 
tation energy was made with reasonable agreement 
with available measured fragmentation cross sections 
(ref. 15). Most of the data to which the model 
was compared were of relatively high energies of 0.9 
to 2.1 GeV/amu. Indeed, the cross sections for frag- 
menting a projectile of mass A p into a fragment A p 
with charge Zp showed relatively little dependence 
on energy (especially for light-fragment production 
(n, p , d , t, he, and a)) as seen in figure 1 for H 2 O 
targets. (Note that the early versions did not cal- 
culate d, t, and he cross sections, which were added 
only recently.) This absence is especially disturbing 
because one would assume that the removal of 80 per- 
cent of the mass from an iron nucleus might require 
more energy than the kinetic energy of 25 MeV/amu. 
Obviously, three important factors are missing from 
the NUCFRAG and HZEFRGl (a publicly released 
version of NUCFRAG) codes. First, the straight-line 
trajectories throughout the interaction assumed in 
NUCFRAG underestimate the closeness of approach 



between the projectile and target. Second, the ki- 
netic energy of the projectile is downshifted as en- 
ergy is given up to release nucleons in the collision 
event. Third, energy conservation occurs in particle 
knockout processes, and this effect further modifies 
the trajectory. The first two factors will be addressed 
herein. 


Theory 

The essence of the NUCFRAG (and HZEFRGl) 
code is that the mass removal A A from a projectile 
of mass Ap by collision with a target of mass Ap 
is a function of the impact parameter b. The as- 
sumption of straight-line trajectories also makes the 
impact parameter the distance of closest approach. 
Swiatecki, Bowman, and Tsang (ref. 11) suggested 
that the cross section for removal of A A nucleons 
er(AA) is given by 


r 

r f i\ j 

a(AA) = 7r ^ 

r ( 2 )] 



which has brought these simple concepts to useful- 
ness. We maintain this relationship but no longer 
take b to be the distance of closest approach. 


Coulomb Trajectories 

The equations of motion in the nuclear coulomb 
field are given by energy conservation as 


where 


r m — 


ZpZpe 2 

Etot 


(6) 


Note that r m is the distance of closest approach for 
zero impact parameter. We now take the distance 
of closest approach r calculated by NUCFRAG at a 
given A A and calculate the impact parameter. Thus, 
we extrapolate backward in time along the coulomb 
trajectory to the initial impact parameter b using 
equation (5). This calculated value of b is used in 
equation (1) to evaluate the cross section. Note that 
the effect of the coulomb trajectory is to move the 
separation at impact r to smaller impact parameters 
b and thus reduce the cross sections, especially at low 
energy. 


Energy Downshift 

A second correction to the trajectory calculation 
comes from the transfer of kinetic energy into binding 
energy in the release of particles from the projectile. 
(Obviously, energy is also lost in releasing particles 
from the target, which we do not yet calculate.) The 
total kinetic energy in passing through the reaction 
zone is reduced to 


E f = Ei- 10 AA (7) 

by assuming that 10 MeV is the average binding en- 
ergy. The kinetic energy used in the closest approach 
calculation is the average 


Em = \t*r 2 + 


2 1 ” ' 2 iir 2 


+ 


ZpZpe 1 


( 2 ) 


Em = \(Ei + Ef) — Ei — 1(10 A A) (8) 


where Etot is the total kinetic energy in the center of 
mass system at large relative distances, r is the rel- 
ative distance between the charge centers with time 
derivative r, p is the reduced mass, l is the angular 
momentum, Zp and Zp are the atomic numbers of 
the projectile and target nucleus, respectively, and e 
is the electric charge. (That is, e 2 — 2R y a 0i where R y 
is the Rydberg constant and a 0 is the Bohr radius.) 
The angular momentum is given as 

e 2 = 2 fiE M b 2 (3) 


Obviously, Etot as given by equation (8) is very crude 
and substantial improvements can be made. Isobaric 
fragmentation cross sections of the old NUCFRAG 
model are shown in figure 1, and the coulomb revised 
cross sections are shown in figure 2. (Note that the 
curves labeled as mass 0-5 are for the light fragments 
n, p, d, t 1 he , and a.) The aforementioned formalism 
greatly improves the nuclear data base, as can be 
seen by comparing figure 1(a) with figure 2(a), but it 
is relatively ineffective at the higher energies above 
600 MeV/amu. (Compare fig. 1(d) with fig. 2(d).) 


The distance of closest approach is given by equa- 
tion (2) for f — 0 as 


__ E to tb 2 ZpZpe 2 
r 2 r 


(4) 


which we write as 

b 2 = r(r - r m ) (5) 


Comparison With Experiments 

As a test of the formalism, in figures 3 and 4 we 
compare our theory with the experimental data ob- 
tained by references 16 and 17, respectively. Quali- 
tative agreement is achieved for charge removal from 
projectiles of O and Kr, even at very low energies. 
This agreement can likely be improved if the inter- 
ference of nuclear dissociation and coulomb phase 
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shifts is calculated. Clearly, transfer processes are 
important at these energies and are not represented 
in the current model. The large cross section for 
AZ = 0 observed for the 16 0 experiment in figure 3 
is probably due to exchange contributions at these 
low energies. 

The elemental cross sections for ie O + 27 A1 re- 
actions are shown in figure 5 at low energies. The 
HZE fragmentation cross sections in figure 5(a) ex- 
perience the greatest variation below 25 MeV/amu, 
an effect arising mainly by coulomb deflection. A 
somewhat weaker variation above 25 MeV/amu re- 
sults from the varying mean free path of the two- 
body collision. The light- fragment-production cross 
sections show rather strong energy dependence below 
200 MeV/amu, as seen in figure 5(b). The elemental 
cross sections for 16 O + 204 Ti are shown in figure 6. 
The much stronger energy dependence shown in com- 
parison with that of the aluminum target in figure 5 
is due to the much greater charge of the titanium 
target. 

Increasing the projectile charge produces an even 
greater energy variation of the HZE fragmentation 
cross sections, as shown in figure 7. The HZE frag- 
mentation cross sections show even greater energy 
dependence, especially for AZ — 0, 1 and AZ > 20. 
The nearly linear dependence of the AZ = 0, 1 cross 
sections is due to the coulomb excitation of the gi- 
ant dipole resonance and the resultant coulomb dis- 
sociation (ref. 18). The strong variations for AZ > 20 
are due to the coulomb distortion of the collision 
trajectory. In all cases, the effects are most easily 
observed in the light-fragment-production cross sec- 
tions which are sensitive to small changes in frag- 
ment distributions caused by the large multiplicities 
of light fragments, especially for heavy projectiles. 

Concluding Remarks 

The nuclear fragmentation cross sections are im- 
proved by calculating the interaction along the clas- 
sical coulomb trajectory. An energy-loss term that 
is proportional to the mass removed in fragmenting 
is used as a correction to the distance of closest ap- 
proach. The results are physically more meaningful 
at low energy than the prior code NUCFRAG (and 
a publicly released version called HZEFRG1), but 
they approach the older results above 600 MeV/amu. 
Further improvements in estimating the effects of en- 
ergy loss should improve the present approximation. 
The adequacy of the present cross sections at low en- 
ergy must await further experimental comparisons. 


Reasonable agreement is obtained with low-energy 

charge-removal cross sections. 

NASA Langley Research Center 

Hampton, VA 23681-0001 
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Figure 3. Charge-removal cross sections for 11.7 MeV/amu 16 0 projectiles onto 92 Mo targets are shown with 
experimental values obtained by reference 16. 



Figure 4. Charge-removal cross sections for 84 Kr projectiles onto 166 Er targets are shown with experimental 
values obtained by reference 17. 
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(a) HZE fragmentation. (b) Light fragment production. 

Figure 5. Elemental cross sections for 16 0 fragmentation on an aluminum target. 
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